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Abstract 

We show that theories with ghosts may have a long hved vacuum state even if all interactions 
are Lorentz preserving. In space-time dimension D > 2, this is only possible if we assume that 
the interaction between ordinary matter and the ghost sector is soft in the UV, so that it can be 
described in terms of non-local form factors rather than point-like vertices. 
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I. INTRODUCTION 



Ghosts have the reputation of being disastrous for particle theories. Any theory that 
allows ghosts, that is, particles with negative kinetic energy, to appear in physical states 
suffers from an instability with respect to spontaneous creation of ghosts and ordinary 
particles with zero total energy and momentum. This process may be strongly suppressed 
if the coupling between the particle and ghost sectors is very weak. But even in the absence 
of any direct coupling both sectors must be coupled to gravity, and so gravity-mediated 
vacuum decay cannot be avoided. 

Furthermore, the following standard argument [l-3| indicates that the corresponding 
decay rate is divergent. Each particle-ghost nucleation event can be characterized by the 
4- momentum P'^ of the created particles, the 4- momentum of ghosts being K'^ = —P^. The 
vacuum decay rate per unit spacetime volume can depend on only through the invariant 

s = P^P^. (1) 

Hence, the total rate will be given by an integral over s and over all possible locations of 
the 4- vector on the mass shell ([1]), 

where F{s) is a mo del- dependent function. The P-integration is over the Lorentz group. 
This integral is manifestly divergent, which indicates that the theory is inconsistent. 

On the other hand, ghosts have been proposed for some interesting cosmological roles. 
Linde [4] has argued that the smallness of the cosmological constant can be attributed to a 
ghost sector described by the same Lagrangian as the matter sector, but with opposite sign. 
Kaplan & Sundrum suggested [sl that the existence of such a sector could follow from an 
approximate discrete symmetry they called "energy parity". This symmetry is broken by 
gravity, and a small cosmological constant is induced. Theories with ghosts have also been 



suggested as models of dark energy jl, |2|, |5|. These ideas provide motivation not to dismiss 
ghosts too lightly and to find a way to get around the vacuum decay catastrophe. 

A possible, if somewhat radical, way to do that was pointed out in Refs. As- 
suming that the interaction between the particle and ghost sectors is mediated by gravity, 
they suggested that both integrations in (|2]) are effectively cut off due to modifications in 
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gravitational physics above some energy scale fi. In order to induce a cutoff in the Lorentz 
group integral over P, this new physics necessarily has to violate Lorentz invariance. With 
these assumptions, the vacuum decay rate can be estimated as 

p 

and it was argued in 21] that for fi MeV particle production at this rate out of vacuum 
is consistent with observations. The expected magnitude of the cosmological constant in 
this kind of models is 

A ~ /i^ (4) 

which agrees with the observed value for fi ~ 10^"^ eV. 

Giving up Lorentz invariance is a high price to pay, and in the present paper we shall 
explore the possibility of living with ghosts in Lorentz invariant theories. We shall argue 
that this is indeed possible. The key point is that even if Lorentz invariance is not explicitly 
broken, a vacuum with ghosts can at best be metastable, so it could not have existed 
forever. If this vacuum was created on some spacelike hypersurface S, this surface breaks 
Lorentz invariance spontaneously and can in principle introduce an effective cutoff for the 
P-integration in ([2]). In a cosmological context, the ghost-carrying vacuum could exist inside 
certain kind of bubbles in models of eternal inflation. In this case, the role of S would be 
played by the hypersurface where the phase transition occurs, from an earlier ghost-free 
vacuum (or, if all vacua are ghost-carrying, then S might be taken as the semiclassical 
beginning of the universe). For simplicity, here we shall restrict attention to a Minkowski 
vacuum, and we shall assume that E is a fiat hypersurface. We shall see that in the presence 
of such a surface the boost integration in (jj]) is indeed regulated, provided that the interaction 
between the particle and ghost sectors which is responsible for the vacuum decay is non-local 
and has a finite characteristic time (and length) scale fi~^. 

Generic forms of on-locality can be problematic in quantum field theories, as they may 
lead to violations of causality. Such violations, however, may be phenomenologically accept- 
able, as long as they are confined to sufficiently small scales. This appears to be the case in 
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Lee- Wick theories, which are nonlocal generalizations of the Standard Model 
problem is that the time ordered exponential in the standard perturbative expansion of the 
S-matrix is not well defined in a nonlocal theory. Time ordering becomes frame-dependent, 
and as a result the Lorentz invariance of the S-matrix is not guaranteed. Difficulties with the 



defininition of a meaningful S-matrix for generic non-local interactions have been emphasized 



m 



a. 



On the other hand, non-locality is acceptable in effective field theories, where it arises 
from integrating out degrees of freedom of an underlying causal theory. In particular, string 
theory is expected to be described by an effective field theory, where non-locality becomes 
important at energies comparable to the string scale. ^ We note also that even though 
a generic nonlocal theory may violate Lorentz invariance, there are some special theories 
which are known to be free of this problem; see, e.g., Refs. {3, liol. We shall assume that the 
nonlocal interaction between particles and ghosts belongs to this special class. 

Our basic idea can be illustrated by the following heuristic argument. Suppose for sim- 
plicity that the ghost-infested vacuum is locally Minkowski and was created on a flat hyper- 
surface t = 0. Suppose further that in the frame where both particles and ghosts nucleate in 
their respective center of mass frames, that is, where P = 0, the spacetime region affected 
by the nucleation process is At ~ |Ax| ~ A nonzero value of P corresponds to the 

same process boosted to a frame moving with Lorentz factor 7 ~ |P|/yU- For 7^1, the 
affected region has the same spacetime volume, but is highly stretched in the direction of 
motion and in the time direction, e.g., 

At~7/i"^ (5) 

For vacuum decay occurring at time t, we should require At ^ t, since otherwise the nu- 
cleation region will stretch beyond the initial surface t = 0. The corresponding momentum 
cutoff is 

Pmax ~ fJ^'^t. (6) 

The particle creation rate is then 

r~^(.*)=. (7) 

p 

Note that integration over the Lorentz group in Eq. ([2]) comes from integration over 
external momenta and does not arise in loop diagrams. Hence, vacuum energy diagrams are 
cut off at /i, not at Pmax, and the gravitational contribution to the cosmological constant is 
still expected to be ~ /x^. 



^ Independent arguments for an intrinsic non-locality of gravity have been presented in the context of black 
hole evaporation and holography. See e.g. and references therein. 
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The possibility that the divergent vacuum decay rate may be regularized by a finite age of 
the universe has been briefiy considered by Cline et al [2]. They suggested that the resulting 
decay rate would be given by an expression similar to Eq. ([7]), where fi has the meaning 
of the high-energy cutoff of the theory. The main difference with our results is that the 
authors of ^ did not recognize the essential role of non-locality in regulating the vacuum 
decay catastrophe. 

In the rest of this paper, we will study the decay rate of a vacuum with ghosts by a direct 
computation, showing agreement with the heuristic estimate ([7]). The paper is organized 
as follows. In Section II we discuss the decay of a ghost-infested vacuum assuming a local 
interaction between ordinary particles and ghosts in D spacetime dimensions. We find that 
a finite vacuum lifetime does regulate the decay rate in D = 2, but not in any higher number 
of dimensions. 

In Section III, focussing now on the physically interesting case of D = 4, we consider 
a simplified model with restricted non-locality. In this model, interactions are local within 
each sector while the interaction mediating between the sectors is described by a non-local 
form factor. This setup is sufficient to eliminate the divergence due to boost integration, in 
agreement with intuitive expectations. However, in this simplified scenario, the form factor 
depends on a single invariant, and this will not be enough to control additional phase space 
divergences which may arise when the number decay products is bigger than two. More 
precisely, we may be left with a divergence due to integration over the invariant energies of 
the decay products within each sector. As we shall see in Section IV, this can be remedied 
by considering fully non-local interactions. In this case, the rate of decay is finite provided 
that the interaction is soft above some energy scale /i, so that all kinematical invariants 
constructed out of the decay products are effectively bounded by that scale. Again, the rate 
agrees with the estimate ([7]). We conclude in Section V. 

II. LOCAL INTERACTION 

Consider a theory with ordinary matter fields 0, and ghost fields ip in D spacetime 
dimensions. A generic vertex with n ordinary particles and n' ghosts has the form 




(8) 
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For simplicity we assume all fields to be scalars. 

To investigate the effect of a finite lifetime of the vacuum, we shall confine the duration 
of the interaction to a finite time interval, 

|x°|<T. (9) 

This can be implemented by means of a Gaussian regulator, 

s.^xj,''..-'i'-r(.)r-(.). (10) 

The interval At ~ 2T will then play the role of the time elapsed since the creation of the 
ghost-carrying vacuum, with the hypersurface t ~ —T playing the role of the surface S on 
which that vacuum was created. (Note that boundaries between spacetime regions occupied 
by different vacua are not sharply defined, so the sign "~" is appropriate here.) 

From ([H]), the matrix element for production of n ordinary particles of momentum and 
m ghosts of momentum ki is given by 

M{0^ Pi...Pn, h, fc„0 = xjd''x e-*(^+^)^ e-^"/^\ (11) 

where P = Y^pi and K = ^ki are the total momenta of particles and ghosts. Integrating 
over d^x, we obtain 

M = {2n)^\ 5^-\P + K) Srin) , (12) 
where Q = + and we have introduced the function 

= e-^^^V4. (13) 

The latter approaches the standard 5{VL) of energy conservation in the limit T — )■ oo. 
The total probability of decay is given by 

V = I n^Pi^\ \M\\ (14) 

where we use the notation 

Note that Ai only depends on the total momenta P and K for particles and ghosts. Hence, 
it is convenient to factor out the phase space integrals over the individual momenta of the 
decay products. To do this, we write 

p = J d^'p d^'K ds ds' 6{P^ - s) S{K^ - s') p^P^) p4K^) \M\\ (16) 



where we have introduced 

/n „ n' 

n rfp. 5^^^ {p - Y.p) ' p^^^') ^ / n 5^^^ - E ^) ■ (17) 
i=i j=i 

From ( !T2|) and f|T6l) . we find the general expression for the inclusive decay rate per unit 
volume: 

^totai = = / dsds' p^{s)p^{s') r(s,s'), (18) 

where 

ns. s) = \{2.rv^ X' j (p2+,)C(p'^+,.)i/2 ^HV2 n) (19) 

and the parameter VL is now expressed in terms of s,s' and P, 

1] = P° + = (P2 + sY'^ - + s'fl\ (20) 

Note that is negative because if) fields are ghosts. 

The rate r(s, s') has to be integrated with the spectral densities p{s) in order to obtain 
the total rate ( fTSj) . For the special case of a single stable particle of mass, n = 1, the density 
is 

p{s) = 5{m'-s), (21) 

where m is the particle's mass. For n > 1, we will be mostly interested in the form of p{s) 
at large values of s, much greater than all particle masses squared. On dimensional grounds, 
we can write 

D(n-l) 

pAs) ~ (22) 

and similarly for Pip{s). 

To connect with standard arguments, let us for a moment consider the limit T — )■ oo. 
Then, 

5t{V2 ^ 5{V2 n) = V2 (P2 + sf'H{s - s'). (23) 
Using f|T9l) and f|T8|) . we have 

Ttotai = ^(2vr)^+5 j ds p^{s)p^{s) r(s), (24) 

where 
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The momentum integral in ( 125|) leads immediately to the well known divergence associated 
with boost integration, as in ([2]). 

Returning now to finite values of T, let us examine the momentum integral r(s, s') in (IT^ . 
At large values of |P|, ^ s,s', we have Q ^ {s — s')/2|P| ^ 1 and 6t{V^^) ~ const. 
Hence, for D > 2 the integral is still divergent. We shall discuss what can be done about 
this divergence in the following sections. But before getting to that, let us stop to consider 
the special case of D = 2, when the integral f lT9|) is convergent. 

In this case, r(s, s') can be estimated as 

T{s,s') ~ X^mm{\s - s'\-^ ,T/V^T^} ~ X^s - s'\ + VTT^ /T)-\ (26) 

To estimate the total rate ( !T8|) . we note that for D = 2 it follows from ( 122|) that the spectral 
densities are p^{s) ~ Ptp{s) ~ for all values of n, n' . Substituting this and ( l26i) in Eq.( fT8|l 
for the total rate, we obtain 

/"°° r/s r°° r/s' , 

Ttotai - AM - / ^(|, - s'l + ^^/TY' ~ X'm-^\n{mT). (27) 

Here, the lower bounds on s and s' integrations are set by the particle masses, which we 
assume for simplicity to have the same order of magnitude, ~ m. Note that the total 
vacuum decay rate f l27|) increases logarithmically with time T, as one might expect from the 
logarithmic divergence of the boost integral ( 125|) in 2D. 

In the case of physical interest, D = 4, the vacuum decay rate f|T9l) is divergent even at 
finite values of T. In the following sections we will show that this divergence can be removed 
if we allow the interaction between the particle and ghost sectors to be nonlocal. From now 
on we shall focus on the physically interesting case of D = 4. 



III. RESTRICTED NON-LOCALITY 



We shall first consider a simplified model in which interactions are local within the particle 
and ghost sectors, while the interaction mediating between the sectors is non-local, being 
described by a Lorentz invariant form factor g{z). The corresponding interaction Lagrangian, 
generalizing Eq.®, has the form 

Si = X I (fx I d^z 0"(x + z) g{z) V^"'(x - z) . (28) 
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A standard local interaction is recovered if we take 

g{z)=5^'\z). (29) 
As before, we restrict the duration of the interaction to a finite time interval, 



[X 



° ± z'^f < T\ (30) 



by means of a Gaussian regulator. Then the dependence on x° and z^ factorizes, as 
g-xg/T^g-zg/yz^ Q]2 the other hand, it may be useful to formally allow for two different 
time-scales T' and T in the two factors of the regulator: 

Si = \ j d^x e"""/^" j (fz 0"(x + z) g{z) e-^o/^V™(a; - z) . (31) 

In the limit T' — )■ oo the interaction is translation invariant, and energy is conserved even if 
we keep T finite. In what follows, unless otherwise stated, we shall keep both of them finite. 
The case of our interest will be T ~ T', corresponding to an interaction of finite duration. 

With a nonlocal interaction (13T|) the matrix element (fT2|) for particle and ghost production 
generalizes to 

M = \ j d'^z g{z) e-^o/T^ j ^A^^-iP{x+z)^-iK{.^z) _ (27r)^A53(P+K) 5t'(^^) Wriw^P) 

(32) 

where 

Wt{w, P) = y" d^z g{z)e-''"''e^'^'e-''/^' (33) 

and we have introduced the notation w = (P^ + s)^/^ + (P^ + s'Y^'^. The total decay 
probability is given by Eq. ffTSj) with 

/^3p 
JWTsMp^TW^^^^^' 5t'(v^ n). (34) 

Expressing the form factor g in terms of its Fourier transform G, 

g{z^) = (27r)-^ J dE d\ e'(^-°-'i^) G{E, q), (35) 
and substituting in have 

/oo 
dE 5t{w - E) GiE^^P). (36) 
-oo 



fOO 

Wt 
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In general, we can split the Lorentz invariant function G{E, q) into negative and positive 
frequency parts 

G{E, 2P) = e{E) G+ + e{-E) Gs, (37) 
where, we have introduced the invariant 



S = {E/2Y - Pi (38) 

Note that G^ should vanish for S < 0, otherwise we would have imaginary energies at low 
momenta (corresponding to exponential growth in the integrand in fl55]) at large time-like 
separations). 

Substituting (136|) in (fTSjl . we have 



(T,(t'=± 



(39) 



where 



Here, 



^aa' ^ (^a / ,3p 5t'{V2 5t{w - E„) 5t{w - E'^,) 

^ ^' / (P2 + s)l/2(p2 + 5/)l/2(p2 + 5)l/2(p2 + ^/)l/2- ^^^^ 



E± = ±2(p2 + 5)^/^ E'^ = ±2{V^ + S'f/^. (41) 

Note that the integral over momenta is manifestly convergent. Further, we note that, unless 
(7 = 0"' = +, the 5t{w — -E'er) exponentially switches off the interaction in the limit T — )• oo. 
As an illustrative example consider the form factor 



Gf = J±5(5-so). (42) 

1/2 — 

For T ^ Sq , the Gg contribution drops out and the rate simphfies to 

Vtotal OC A J+T y ds ds P^{S)P^{S ) J d P ^p2^^^i/2(p2^,/)l/2(p2^,^) - (^3) 

where Eq = 2(P^ + so)^^^. The momentum integral converges, and is dominated by momenta 
of order ~ (s + s')T. Hence, we can estimate 

Ttotai ~ X'^llTT' J ds ds' p^{s)p4s'){s + s')-\ (44) 

We thus see that the quadratic divergence in the decay rate which occurs in ([2]), due to 
integration over the total momentum P of the decay products, disappears in the present 
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context. Instead, we obtain a rate that grows quadratically in time, in agreement with the 
heuristic arguments presented in the Introduction. It can be verified that these conclusions 
are not restricted to our specific example fl42l) and apply to a wide class of form factors. 

However, we should also consider the integral over the invariant energies s and s' of the 
decay products. When s and s' are large compared with the masses of the species involved, 
the spectral densities are given by (cf. Eq. (!22|) with D = 4) 

p^{s)^s^^-\ p4s') isT'-' (45) 

Hence, for n + n' > 2 the integrals over s and s' in (14^ would not converge unless we impose 
some kind of high energy cut-off. ^ 

Let us compare the present situation with the case discussed by Kaplan and Sundrum j^, 
where the integral over total momentum P is regularized by means of a Lorentz breaking 
cut-off. They noted that the integral over the invariant energy s of the decay products 
would also be divergent and pointed out that this divergence could be removed by new 
(Lorentz invariant) physics which might soften the interaction at energies above some scale 
So (which they took to be of the same order as the Lorentz breaking cut-off). Here, instead, 
we find that even if the interaction only carries invariant energy of order sq, we still have a 
divergence when we integrate over the invariant energies of the decay products s and s'. This 
may seem surprising, since we are used to thinking that the invariant energy of the mediator 
should be the same as that of the decay products. However, in the present case energy is 
not conserved because the interaction is switched on and off. As a result, at momenta of 
order ^ (s + s')T, the 5t(^) of "energy conservation" allows s and s' to be very different 
from each other, and to be much larger than sq. 

In the following Section we explore a model where the interaction vertex between par- 
ticles and ghosts is fully non-local, so that the form factor depends on more than just one 
kinematical invariant. We shall also assume that interactions become soft when the values 
of these invariants exceed the scale /i. As we shall see, this may lead to acceptable decay 
rates in phenomenologically interesting scenarios. 



^ The rate would be finite without any cut-off in the invariant energies if we assumed that = 0, but 
then the two sectors would be disconnected at T — oo, where driw + Eq) ^ 0- 



11 



IV. FULLY NON-LOCAL INTERACTION 

Here, we consider non-local vertices where the different legs of matter in a given sector 
do not necessarily meet at the same point. For simplicity, we restrict attention to the four 
leg vertex: 

Si^xJ d^x d'^zd^yid^y-2(t){x + z + yi)4){x + z-yi)g{z,yi)i^{x-z + y2)i^{x-z-y2). (46) 

The regulated amplitude is then given by 

M = (27r)^A 6^(P + K) 5t'(^) Wt{w, P,p, k), (47) 

where, as usual, P — Pi + P2, K — ki + = P° + K'^, w — — K^, and we have 

introduced p = Pi — P2, k — ki — k2 and 

WT{w,V,p,k)^ j d'zd'y,d'y2 ^(^, y,)e-(-^°-^P^)e-^(^•^^+'=■^^)e-[(^°)'+(^?)'+(^^°)'l/^^ (48) 

In terms of the form factor in Fourier space G, defined by 

g{z, y,) = (27r)-^^ J d\ d\ d\ e^(^-+«.-^i+«^-s'^)G(g, q„ q,), (49) 



we have 

Wt 



/oo 
dE dEp dEk 5t{w - E) Srip" - E^) 5t(A;° - Ek) G{q, q^, qk), (50) 
-00 



where in the last equation q = {E,2P), qp = {Ep,p) and q/. = {Ek,k). Introducing the 
invariants 

S = q^/A = E^/A - P2, (51) 

S, = = El - p^ (52) 

Su^ ql^ El - k^ (53) 

we may decompose G into negative and positive frequencies for each one of the energies. 

G{q.q,.qu)= J] 0{aE)e{a' Ep)e{a" E,)G''^^^;;^^ (54) 

cr,cr',o-"=± 

Then, we may write 

_ 1 ^ / h{w - E^) Srif - E;') 5T{k^ - Ef) 

= A^/^^J "^^^ "^^^ (p2 + ^)i/2(p2 + ^^)i/2(k2 + ^^)i/2 GZ:,,- (55) 
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Note that G'^"^^'^^ depends on the subscripts only trough invariants constructed out of the 
momenta g, qp and qk- These are S", S'p, Sk and the products 

7i"p = g ■ gp, TTfc = g ■ g^, Upk = qp ■ qk- (56) 
The decay rate is given by 

r ^ 

Ttotai = Svr^v^A^ / J] dp, dkj 5\P + K) 6T'iV2 n) \Wt\^. 

i,7 = l 



(57) 



We have (3 — l)x4 — 3 = 5 powers of momentum from the phase space integrals. On the 
other hand, there are 6 powers of momentum in the denominator of Also, assuming 

that the non-local interaction becomes soft at a short distance scale the domain of the 
integration in fl55|) is effectively finite S,Sp,Sk ^ /x^, and we do not expect divergences in 
the decay rate. 

Let us define = yU^T, and assume that T ^ fi^^. Then, we can estimate 

Wt ^ G{q,qp,qk). (|pi| , |ki| < p,) (58) 

Then, assuming that the rate fl57|) is dominated by momenta of order p^,, we have 

Ttot. - T j^^^^ ^0 ,0 feO| MP + K) |G(g, g„ g,)| . (59) 

As mentioned above, the amplitude G, depends in principle on the 6 invariants we may 
construct with g, g^ and g^. Assuming that the interaction becomes soft at the energy scale 
/i, we may require S, Sp, Sk, vr^, iik, T^pk ~ yU^- The first three of these constraints will restrict 
the values of the energies E, Ep, Ek to be very close to the modulus of the corresponding 
spatial momenta. For instance, 



E = 2VP2 + S = 2|P| + 0(5/|P|) ^ 2|P| + 0(1/T). (60) 

where in the last step we have assumed S ^ fi"^ and |P| ~ = fi^T. Similarly, from 
Sp, Sk ^ /i^ we have 

Ep^\p\ + 0{1/T), Ek^\k\+Oil/T). (61) 
Now, from the condition 7Cp ^ we have 

np = q-qp = E Ep-2Pp^ 2|P||p| - 2Pp + 0(p,/T) < (62) 
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The last constraint implies that P and p have to be aligned within an angle 

dp ^ n/p*. (63) 

A similar argument can be made for Tr^ and iTpk, leading to additional constraints in the 
alignment of spatial momenta 

Ok,dpk ^ fJ'/p*- (64) 

To satisfy (163!) and p and k must remain within a solid angle Afip, Afifc ~ '^f^'^/pl 

of the total momentum P, effectively suppressing the phase space available by a factor of 
TT^/i^/p^ (while the integration over the total momentum P spans the full solid angle Air). 
From this, the rate fl5^ can be estimated as 

r«., ~ x'TV: giGP ~ 10-^1 (^r)^ (65) 

in agreement with the heuristic estimate ([7]). In the last step we assume T ~ T' and an 
interaction of gravitational strength AjGI ~ fi'^Mp'^, where Mp is the Planck mass. 



V. DISCUSSION 



Vacuum is unstable in the presence of ghosts, and formal calculations give a divergent 
rate of vacuum decay. Here we found that the rate may actually be finite: the divergent 
momentum integration may be effectively cut off due to the finite lifetime of the metastable 
vacuum. We found, however, that this mechanism can work only if the interaction between 
the particle and ghost sectors is nonlocal above some characteristic energy scale /i. The 
momentum cutoff Pmax is then 

Pmax(t)~A (66) 

and assuming that the ghost-particle interaction is gravitational in nature, the vacuum decay 
rate is given by Eq. f l65|) . 

^totaiit) ~ i^^tf, (67) 

■p 

in agreement with the heuristic argument given in the Introduction. Here, t is the time 
elapsed since the creation of the vacuum, and A is a numerical phase space factor, which we 
estimated to be of order A ~ 10^^. 

Throughout the paper we assumed for simplicity that the metastable vacuum was created 
on a fiat hypersurface t = in Minkowski space. More realistically, such a vacuum could 
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be formed as a bubble in the course of eternal inflation. The role of the vacuum creation 
surface would then be played by a hyperbolic 3-surface of constant negative curvature. For a 
Minkowski vacuum inside the bubble, we expect that the same estimate f l67p for the vacuum 
decay rate should apply, with t now being the proper time in the FRW coordinates inside 
the bubble. For a de Sitter vacuum with energy density > 0, this estimate should still 
hold for t < H-\ where = {SnGpjSY^'^ is the Hubble expansion rate in that vacuum. 

Is it conceivable that our own vacuum is unstable with respect to production of particles 
and ghosts? In such a scenario, the time t in Eqs. ( 166|) .( 167|) should be identified with the 
present cosmic time, to ~ 10^" yrs, and the characteristic energy of the particles being 
produced at present is given by 

E ~ PmaM ~ 10-V-3^P, (68) 

where /i_3 = fi/10~^ eV. 

The total energy density of particles (and ghosts) produced up to the present time is 

~ r(to)Pmax(to)to ~ A^(^to)'. (69) 

p 

This constitutes a fraction 

^ (!_)'(,,„).. ^(^)^10-M.i^3 (70) 

of the total energy density of the universe, po ~ ^p/^o- 

As we shall see, the cutoff scale p must be rather low. This means that to the lowest order 
in pertrubation theory, only massless (and very light, m ^ p) particles can be produced.^ 
Unstable particles will eventually decay into gamma rays, protons, electrons and neutrinos, 
and possibly into some presently unknown stable particles. High-energy photons will scatter 
on the CMB, producing e+e~ pairs and cascading down in energy. The resulting back- 
ground of cascade photons with energies ^ 100 GeV is constrained jisl . 1^ by Fermi-LAT 



measurements 



15| to have energy density less than 

10"^ eV/cm^ ~ 10-%. (71) 



max 
cas 



Heavy particles, of mass m ^ ii. may be produced by multiple graviton exchange. In particular, semi- 



classical pairs of Bondi dipoles may be produced non-perturbatively The instantons describing 



such processes are Lorentz invariant, and as a result the corresponding decay rates do not contain boost 



divergences: the rates are finite without the need of regulating the duration of interaction (see e.g. [l2[). 
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Assuming that the cascade radiation accounts for a substantial fraction of the total energy 
density fl69|) of the produced particles'^, the corresponding bound on /i is 



/i < (1.8 - 5.6) ■ 10"^ eV, (72) 

where the lower value corresponds to A ~ 1 and the higher one to A ~ 10^^. 

On the other hand, we should make sure that the nonlocal modification of gravity we 
assumed here is consistent with the precision short-distance tests of the Newton's law. As 
of now, the inverse square law for forces of gravitational strength has been tested [l^ down 
to distances / ~ 0.1 mm. This corresponds to the energy ~ 2 • 10~^ eV, and so we should 
require 

/i > 2 • 10~^eV, (73) 

which is marginally consistent with ( l72ll . 

As we mentioned in the Introduction, a nonlocal interaction generally leads to violations 
of causality. With fi ~ 10^^ eV, such violations should occur on relatively large distance 
scales, up to 0.1 mm. For electroweak strength interactions these violations could manifest 



themselves in scattering experiments 



17| . However, for a gravitational strength interaction 



we seem to be in no danger of running into conflict with observational bounds. 
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